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8 Abstract: Cotton (Gossypium hirsutism) fiber, an extremely elongated and thickened single cell 

9 of the seed epidermis, is the world's most important natural and economical textile fiber. Unlike 

10 Jrah'dopsis\Qdii trichomes, fiber initials are randomly developed and frequently form in adjacent 

1 1 seed epidermal cells and follow no apparent pattern. Numerous publications suggested cotton 

12 fiber development shares a similar mechanism with Jrafo'c/opsis\Q2& trichome development. 

13 Here we show that H202 accumulation in cotton ovule epidermal cells by NBT staining ovules 

14 at different development stages between TM1 and Nln2, a lintless-fuzzless doubled mutant 

15 originated from TM1. In contrast, Jrafo'cfopsis 'and cotton leaf trichomes do not show H202 

16 content. By adding DPI (H202 inhibitor) and SHAM (H202 activator) in vitro ovule cultures, 

17 we show fiber initiation directly involves with H202 accumulation. We propose that the 

18 directional accumulation of H202 in cotton ovule epidermal cell is the drive for fiber initiation, 

19 elongation. 

20 
21 
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1 Main Text: 

2 Seed dispersal is sometimes split into autochory (when dispersal is attained using the 

3 plant's own means) and allochory (when obtained through external means). Wind dispersal 

4 (anemochory) is one of the autochories and more primitive means of dispersal and plant seeds 

5 have a feathery tissue attached to their seeds and can be dispersed long distances. Primitively, 

6 plants distributed the seeds majorly through autochory; as animals emerge and human activity, 

7 some plants lost their ability to produce seed secondary tissues helpful for dispersal and some 

8 plants gained the production for specific seeds tissues such as cotton seed trichomes as human 

9 activity. But we do not understand how these tissues develop. Cotton {Gossypium spp.J fibers are 

10 unicellular and single-celled seed trichomes of economic importance (7) and provide a special 

1 1 single-cell functional genomics research platform {2, 3) for fundamental biological processes 

12 such as organogenesis and cell wall development in plants {4, J). Fiber development initiates on 

13 the day of anthesis (0 dpa) and 30-40% of primordial cells (d) on the ovule epidermis form 

14 cuboidal outgrowths called fiber "initials" (/). Unlike Arabidopsis\Qdl trichomes, fiber initials is 

15 randomly distributed and follow no apparent pattern as fiber initials frequently form in adjacent 

16 epidermal cells (S,9). Fiber initials rapidly elongate to produce highly tapered cells within 2-3 

17 days of anthesis through a biased diffuse growth mechanism {JO). 

18 The initiation of single celled cotton cuboidal fiber cells is a biologically important and 

19 complex process, but how the fiber initiate and how the genes involved are regulated remain 

20 elusive. Zhang et al., (//) demonstrated by adding 30% H202 to in vitro cultured ovules of 

21 XinFLM, a linted-fuzzless mutant, that the fiber initiation is partially rescued, but no evident 

22 effect on Nl, the dominant naked seed mutant from TM-1 genetic background. This implies the 

23 H202 biosynthesis pathway may be defect and the signal transduction is still functional in 
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1 XinFLM. Like cuboidal outgrowth fiber initials, Andriunas et al., (J2) shows a flavin-containing 

2 enzyme (NAPDH oxidase) and superoxide dismutase cooperatively generate a regulatory H202 

3 to govern ingrowth wall formation, a transfer cell (/J), using Visia Faba cotyledon culture. 

4 Forman et al., (J4) verified that ROS accumulates in growing wild type root hairs but 

5 dramatically decreases in rhd2 short hair mutants and the generation of ROS is inhibited and 

6 mimics rhd2 mutant and further showed (/J) a RhoGTPase GDP dissociation inhibitor controls 

7 spatial ROS production during root hair development. 

8 To determine whether H202 is necessary for cotton fiber initiation, we firstly used SEM 

9 to observe fiber initiation and development dynamics starting ovules at 2 days of pre-anthesis to 

10 2 days of post-anthesis. In the dominant Nl mutant (Fig. 1), no fiber initials are present on ovule 

1 1 at anthesis and very rare fiber initials occur on the ovules at 2 days of post anthesis. The 

12 recessive mutant n2 does not show significant change for fiber initials. The double mutants Nln2 

13 generated by crossing Nl and n2 and selecting till 7 generation does not show any fiber initiation. 

14 XZ142-fl, a lintless and fuzzless mutant originated from its wild type XZ142, also does not 

15 produce any fibers. 

16 In higher plants ROS can probably be produced by several alternative pathways (/d) 

17 including a cell-wall-localized peroxidase, amine oxidases and non-flavin NADPH oxidases, or 

18 NADPH oxidases that resemble the flavin-containing type that is activated by Rac in leukocytes 

19 (//}. Simultaneously, a vast network of antioxidants is acting like ROS scavengers of H202 and 

20 maintains ROS homeostasis (73). This scavenging system consists of catalase (CAT), ascorbate 

21 (APX) and secretory peroxidases (POX), glutathione reductases (GR) and peroxiredoxines 

22 (PRX), and non-enzymatic compounds like tocopherols, ascorbic acid and flavonoids. (J 9). Li et 

23 al., (20) identified a cotton cytosolic APX1 (G/iJPJTf) to be highly accumulated during cotton 
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1 fiber elongation by proteomic analysis and showed {21) ethylene is induced in in vitro cultured 

2 ovules by endogenous H202 and found a feedback regulatory system between H202 and 

3 ethylene in fiber development. Hovav et al., {22) found three genes modulating hydrogen 

4 peroxide levels were consistently expressed in domesticated and wild cotton species with long 

5 fibers, but not detected by q-PCR in wild species with short fibers. Potika et al., {23) Potikha et 

6 al., 1999) speculated H202 may function as a developmental signal in the differentiation of 

7 secondary walls in cotton fibers with support of coincidence of H202 generation with the onset 

8 of secondary wall deposition and prevention the wall differentiation because of inhibition of 

9 H202 production or removing the available H202. H202 modulates downstream signaling 

10 events including calcium mobilization protein phosphorylation and gene expression {24) which 

1 1 may lead to cotton epidermal cell fate determination. To determine whether H202 involves fiber 

12 initiation, we used O.lmg/mL NBT (Nitro Blue Tetrazolium), a chemical specifically staining 

13 H202 in living cells (Fig.2). Fig.2A timely profiled a normal cotton fiber development starting 5 

14 days of pre-anthesis to fiber maturation at 50 days of post-anthesis. Fig.2B demonstrated that 

15 TM1 ovules stained darkly at 2 days of pre-anthesis to 2 days of post-anthesis has already 

16 accumulated H202 in comparison to these from Nln2. Significantly, in contrast to TM1 ovules, 

17 no any staining was detected on Jrabidopsis and Cotton leaf trichomes (Fig.2C). We also used 

18 cut cotton leaf edge as positive, dark staining was found at the cut leaf edge and spots with high 

19 concentration of polyphenols (Fig.2C). To rule out the possibility of staining accessibility, we 

20 conducted H202 production measurement using an Amplex red hydrogen peroxide/peroxidase 

21 assay kit (Molecular Probes, http://probes.invitrogen.com/). We found significant decrease of 

22 H202 in ovules of different stages of Nln2 comparing with those from TM, especially at the 2 

23 DP A (Fig.2D). 
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1 To verify whether H202 directly influence epidermal cell differentiation and initiation, 

2 we conducted in vitro ovule cultures by adding DPI (H202 inhibitor) and SHAM (H202 

3 activator. The result (Fig.3A) shows the fiber initiation normally developed in the control (with 

4 0.1% DMSO) in vitro ovule culture and NTB staining shows H202 accumulation in the fiber 

5 initiation and epidermal cells. Fig.3B shows the fiber initiation was almost totally retarded in the 

6 in vitro ovule culture with addition of 0.3 uM DPI and the NTB staining shows barely H202 

7 accumulation. Interestingly, Fig.3C shows very similar results with Fig.3A. We cannot conclude 

8 SHAM enhances fiber initiation as a more controlled experiment is needed for this result. Based 

9 on this data, we can efficiently conclude that H202 directly involves cotton fiber initiation and 

10 development. 

1 1 To molecularly characterize the mechanism of fiber initiation, many researchers {25, 2d) 

12 proposed the mechanism similarity of cotton fiber development to Jrafo'dopsisftvcfaomg, and 

13 study cotton fiber initiation by cloning cotton orthologos for Jmfo'cfopsis trie home development. 

14 Three criteria were used to characterize whether the genes of interest relate to cotton fiber 

15 initiation including (i) orthologs to Jrah'dopsisixv^aovaQ related genes; (ii) highly expressed in 

16 cotton ovules and fiber initials and (iii) specific expression in Jrabidopsis transgenic Col-0 

17 trie home by cotton orthologous promoter-driven GUS and complementation of trichome 

18 development for Jrafo'dopsisixi^omsX&ss, mutants by overexpression of cotton ortholog and the 

19 varied influences to cotton fiber development in loss-of-function and gain-of function transgenic 

20 cotton {27). Not like Jrafo'dopsis trichome, cotton fibers are initiated from seed epidermal cells, 

21 their initiation and development highly synchronize and unbranch with extreme elongation {2S). 

22 To verify the reasonability of the speculation, we roughly grouped the documented cotton fiber 
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1 development related genes into three groups (S-table 1) including ROS-related, fiber initiation- 

2 related and fiber elongation-related. 

3 Transcript profiling for ROS-related genes including G/im/rSPS, with its two targets, 

4 GACSDSmd G6CSD2, GAAPJT, GhRDLJwA GMDL2, only GhRDlIwA GABZ?Z2show 

5 consistently lowest expression in the two lintless fuzzless mutants but significantly increase 

6 expression starting 2 days of pre-anthisis and post-anthesis (fig. IS). The expression of these 

7 genes helps maintain homoestatisis of H202, an important signal for cotton fiber initiation and 

8 development. We grouped cotton orthologs based on homology to Arabidopsis leaf trichome 

9 development into cotton fiber-initiation-related genes. We do not find specific correlation of 

10 these group genes to the two lintless fuzzless mutants. Interestingly, all of these genes express in 

11 the cotton expanded leaves which produce trichomes (fig.2S). This result suggests these genes 

12 may be cotton leaf trichome orthologs but not cotton seed fiber genes. 

13 Interestingly, Walford et al., report {29) RNAi cotton transgenic plants for G/iMFB2J-/ike 

14 resulted in fibreless seeds but developed with normal trichomes elsewhere, mimics the Xul42 fl 

15 mutant. Furthermore, these plants had decreased expression of the fibers related MYBs including 

16 GAMFB2S and G/iMFBJOP, suggesting these MFBs are downstream of G/iMFB2S-Me. 

17 G/iMFB2J-/ike, not an ortholog to Arabidopsis trichome genes, is identified from its reduced 

18 expression in a fibreless mutant of cotton (Xul42 fl) (J^Du et al., 2001). These data implicate 

19 cotton fiber initiation and development may share similar regulatory process with Jrafo'cfopsis 

20 trichome but the involved genes are different. It is important to note that Arabidopsis seeds do 

21 not have seed trichomes like cotton and thus cannot be considered as an ideal surrogate model to 

22 evaluate the specificity of expression of cotton fiber genes. The search for cotton fiber initiation 

23 and development genes is still elusive. We also detect the expression of the third group of 
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1 documented fiber related gene, fiber elongation genes, lowest to undetectable expression of 

2 G/iMFB2J G/iBAPJ (fig.3S) relate to the two lintless fuzzless mutants. This implies these 

3 two genes are acting downstream and necessary for fiber elongation. 

4 Cell morphogenesis is the process by which cells acquire shape during development and 

5 is central to organogenesis in multicellular organisms. Based on the currently available results, 

6 Qin et.al. (37) speculated that fiber cells may elongate via a combination of both tip-growth and 

7 diffuse-growth modes and named linear cell-growth mode. Using NBT staining ovules at 

8 different development stages between TM1 and Nln2, a lintless-fuzzless doubled mutant 

9 originated from TM, we show that H202 accumulation in cotton ovule epidermal cells but no 

10 H202 staining in Jrab/c/ops/s 'and cotton leaf trichomes. We propose that the directional 

1 1 accumulation of H202 in cotton ovule epidermal cell is the drive for cotton fiber initiation and 

12 elongation. 

13 
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1 Fig. 1. Altered fiber initiation and development in two pair of naked seed mutants with their wild 

2 types. (A) Almost all fiber initials occur at 0 DPA in wild type TM1. (B) Few fiber initials occur 

3 at 0 and 2 DPA in dominant mutant Nl . (C) Many fiber initials occur at 0 and 2 DPA in receive 

4 mutant n2. (D) No fiber initials developed on JVJn2 double mutant. (E) Almost all fiber initials 

5 occur at 0 DPA in wild type XZ142. (F) No fiber initials produced on ovules of XZ142-A. Black 

6 Bar = 200um. 
7 

8 

9 Fig.2 Alteration of H202 accumulation by NBT staining and measurement. (A) The time scale 

10 indicates the duration of epidermal cell initiation, fiber initial elongation, cellulose biosynthesis 

11 and maturation. (B) Ovules from TM1 and Nln2, stained with NBT (C). Leaf trichomes from 

12 Arabidopsis and TM1 and TM1 leaf cut edge, stained with NBT. (D). More H202 accumulation 

13 in ovules from TM1 at different stages comparing to those from Nln2.Red Bar=lmm. 
14 

15 Fig.3. In Vitro Ovule Cultures Treated with DPI and SHAM and Stained with NBT. 

16 Cotton ovules from TM1 plants were collected at 0 DPA, treated with DPI (inhibitor at 0.3uM) 

17 or SHAM (activator at 1.5mM) dissolved in DMSO separately for 48 Hrs, and stained with NBT 

18 for 2.5 Hrs. Red Bar=lmm. 
19 
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1 Materials and Methods 

2 Plant Materials and Scanning Electron Microscopy 

3 Upland cotton (Gossypium liirsutum Z. cv TM-1), its isogenic naked seed mutants (N1N1 and 

4 n2n2), the double mutant Nln2 and XZ142 with it fiberless fuzzless mutant XZ142-A were 

5 grown in the greenhouse under 32/2 1°C (day/night) temperature regime. Ovules were collected 

6 at -2, 0 and 2 dpa flowers at 9 and 10 A.M. Hitachi tabletop TM-1 000 scanning electron 

7 microscope was used for imaging the fresh ovules. 

8 Grouping and OF T-PCF Primers ofFotential Genes Felated 'to Fi6er Development 

9 We grouped the potential genes related to fiber development into three categories 

10 including ROS-related, fiber initiation and elongation (Table-Si). 

11 Imaging and FOS Staining for 7M/andJVfn2 Ovules at -2D PA to 2DFA 

12 Ovules from TM1 and Nln2 at -2, -1, 0, 1 and 2 DPA were covered with freshly made 

13 NBT (nitroblue tetrazolium) solution (in 0. 1M potassium phosphate at pH7). Plates were left for 

14 2.5 firs at room temperature for the color to develop. Ovules were imaged under bright-field 

15 illumination using an Olympus SZH10 Zoom Stereo Microscope. 

16 J/2Q2 Measurement 

17 To measure H202 concentration, ovules at different stages were ground in liquid 

18 nitrogen, and 200 mg of ground tissues from each sample was added to 200 uL of 20 mM 

19 sodium phosphate buffer (pH 6.5) and mixed well and centrifuged at 9500g for 10 min. at 4° C. 

20 The supernatants were used for assay. H202 concentration was measured using an Amplex red 

21 hydrogen peroxide/peroxidase assay kit (Molecular Probes: http://probes.invitrogen.com/). 
22 

23 
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1 In Fttro Ovu/e Culture and Treatment m'tn H202 In/ii&otr and Activator 

2 Cotton ovules from TM1 plants were collected at 0 DPA, soaked in 70% ethanol for 1 

3 min, rinsed in distilled and deionized water and placed in liquid media formulated by Beasley 

4 and Ting {S3) in a 6-well plate under aseptic conditions. The ovules were cultured in vitro, as 

5 described by Meinert and Delmer {SP). DPI (inhibitor at 0.3uM) or SHAM (activator at 1.5mM) 

6 dissolved in DMSO was added to medium separately. The final concentration of DMSO was 

7 0.1% in all of the samples. 

8 JUVA Preparation artd 0HT-PCR 

9 Total RNA was extracted from ovules, young expanded leaves using the Spectrum™ 

10 Plant Total RNA Kit and on-column DNase I digestion, following the manufacturer's 

11 recommendations (Sigma- Aldrich). RNA was quantified using nondrop- 1000 and the integrity of 

12 total RNA was confirmed by denaturing agarose gel electrophoresis. QRT-PCR amplification 

13 reactions were performed using SYBR Green detection chemistry and run in triplicate on 96- 

14 wells plates with the 7500 Real Time PCR system (Applied Biosystems). Reactions were 

15 prepared in a total volume of 20 |j.L containing: 5 jxL of template, 1 \iL of each amplification 

16 primer (lOuM), 10 of 2* Power SYBR Green PCR Master Mix (Applied Biosystems) and 3 

17 jxL PCR-quality ddH20. Cotton Histone-3 (AF024716) was used to normalize the amount of 

18 gene-specific RT-PCR products. All PCR reactions were conducted in three replicates and 

19 monitored the absence of primer dimers using a dissociation analysis. The amplification data 

20 were analyzed using ABI7500 SDS software (version 1.2.2), and the relative expression levels to 

21 expression of Histone-3 (fold changes) were calculated. For detecting expression of 

22 GAmiJWASPS, protocol followed (J7) using QRT-PCR primer GTCGTATCCA 

23 GTGCAGGGTCCGAGGTATTCGCACTGGATACGACAAGGGG and forward primer 
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1 ACATTCAAGTGTGTTCTCAGGTCA and universal QRT-PCR miRNA reverse primer 

2 GTGCAGGGTCCGAGGT. 

3 Fig.Sl. Transcripts profiling for ROS homeostasis-related genes. (A) Gene expression for 

4 G/imir398a+b'm different varieties and development stages. (B) Gene expression for 

5 G/iCSDJ(S2J 'm different varieties and development stages. (C) Gene expression for G/iGSD2 \S2) 

6 in different varieties and development stages. (D) Gene expression for G/iJPJfJ '(J 6) in different 

7 varieties and development stages. (E) Gene expression for G/iRDLJ (26) in different varieties 

8 and development stages. (F) Gene expression for G/iJlDL2 '(2d) in different varieties and 

9 development stages. 

10 Fig.2S. Transcripts profiling for potential fiber initiation related genes. (A) Gene expression for 

1 1 G/iMFB2 '(2d) in different varieties and development stages. (B) Gene expression for GfiGPG(S3) 

12 in different varieties and development stages. (C) Gene expression for G/iHOJfJ(2dJ'm different 

13 varieties and development stages. (D) Gene expression for G//J70JT2 '(2d) in different varieties 

14 and development stages. (E) Gene expression for GAHOJO ' (2d])'m different varieties and 

15 development stages. (F) Gene expression for G/iDFL6J(S4) in different varieties and 

16 development stages. 

17 Fig.3S. Transcripts profiling for potential fiber elongation related genes. (A) Gene expression for 

18 G/iMKB2J(29) in different varieties and development stages. (B) Gene expression for 

19 G/iMZB2J-L (29) in different varieties and development stages. (C) Gene expression for 

20 G/iPDFJ(S3) in different varieties and development stages. (D) Gene expression for G/iMZBJ09 

21 fStfJ'm different varieties and development stages. (E) Gene expression for G/iFJTPJ (ST) in 

22 different varieties and development stages. 

23 Table SI. QRT-PCR primers for potential fiber initiation and development genes 
24 
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1 

2 Table SI. QRT-PCR primers for potential fiber initiation and development genes 
3 



Categories 
ROS-related 



initiation 



elongation 



4 
5 
6 
7 



Genes 


NCBI-ID 


Forward primers 


Reverse primers 


G/iCSZ)/ 


F4 15203 


TCTTCTCCCAAGAAGGAGATGG 


TTTGTCGGTGATGGAGAAGCTA 


GACSD2 


EU597268 


TGCGGGTGACCTGGGAAACATAAT 


ATGCTCTTCCAACAACCGCATTGG 


GAJFJri 


EF432582 


TGATTGGTCGACTCAAAACG 


GCCCCAATATTTCCCATCTT 


GAJI/J22-IJ 


AY072821 


AACAAACCCCAATGCAAAAG 


TTGGGTTCCATGCTGATGTA 


G/iJtD22-£2 


AY641991 


CAGCTGGAGTGCAGAAGATG 


GTGTGACAGACTGCTACTGC 


GaMFS2 


AY626160 


ATGGCTCCAAAGAAGGATGGAGT 


GCCGCATCGATTCAAACCTGA 


G/iGPC 


FJ402930 


TGGCTGACTCTCAACATTCCT 


GTTCTACCGGGGATTCTTCC 


GaJ/OJTJ 


AF530913 


GAAATTGGCCCAGAGAATGA 


GAGTCGGAGAGACAGGCAAC 


GaffOJT2 


AF530914 


TGCGAAAGGGTTGCTAGTCT 


ATTGGGTTGACCTGCTTCAC 


GaJ/OJTS 


AY626159 


CCTTCAAAGGATGTGCGAAT 


ATTCGAACCGGAAAGTGTTG 


GAD££6S 


AF336280 


CGAAAGGTCGTTCCTCTTTG 


GCTTCAAGCCATTTTTCTGC 


G/iPDFJ 


AF141375 


GTCCCCCAACTACTCCCATT 


CAAGCTCATGCTTGTTCCAA 


GAMFS2J 


AF336283 


CAATGGGAGAGTGCTCGTTT 


CAACGTTGACGTAGGGGACT 


GAMFB2S-1 


HM134081 


GAGAAATCGAGCCAAGTTGC 


ATCCCCAGAATCACAAACCA 


G/iMrBJ09 


AY072821 


TGGAACAACCACTGAAACCA 


TGGTCAATCATGGCAAGTGT 


GMXPJ 


AF043284 


TCCAGAACAGGATGGCTACC 


AGTTTGTCCGAATTGCCAAC 



8 



